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1 Introduction 

A basket option is an exotic option whose payoff depends on the value of a portfolio of assets. 
Basket options are in general difficult to price and hedge due to the lack of analytic character- 
ization of the distribution of the sum of correlated random variables. Monte Carlo simulation 
is often used to price basket options, which is simple, accurate, but time-consuming. There has 
been some extensive research recently for fast and accurate pricing methods. 

Most work in the literature assumes that underlying asset prices follow geometric Brownian 
motions. The basket value is then the sum of correlated lognormal variables. The main idea 
of the analytic approximation method is to find a simple random variable to approximate the 
basket value and then to use it to get a closed form pricing formula. The approximate random 
variable is required to match some moments of the basket value. Levy (1992) uses a lognormal 
variable to approximate the basket value with the matched first and second moments. The 
results are remarkably good but there is no error estimation. Curran (1994) introduces the idea 
of conditioning variable and conditional moment matching. The option price is decomposed into 
two parts: one can be calculated exactly and the other approximately by conditional moment 
matching method. Rogers and Shi (1995) derive the lower and upper bounds. Vanmaele et 
al. (2004) suggest a moment matching comonotonic approximation for basket options. See Lord 
(2006) for other methods and references. Efforts have been made to extend to more general asset 
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price models. Albrecher and Predota (2004) discuss the NIG Levy process, and Flamouris and 
Giamouridis (2007) the Bernoulli jump diffusion model. 

Xu and Zheng (2009) suggest a jump-diffusion model for underlying asset price processes. 
The innovative feature of that model is that, apart from correlated Brownian motions, there are 
two types of Poisson jumps: a systematic jump that affects all asset prices and idiosyncratic 
jumps that only affect specific asset prices. Such a model can characterize well the market-wide 
phenomenon and individual events. Xu and Zheng (2009) use the partial exact approximation 
(PEA) method to find a closed form approximate solution which is guaranteed to lie between 
the lower and upper bounds. The numerical tests show that the PEA method has superior 
performance in comparison with other methods such as the lower bound, reciprocal gamma, and 
lognormal approximations. The limitation of the PEA method is that it depends crucially on 
the conditioning variable which is derived from the estimation of the basket value and the closed 
form solutions of individual asset prices. This may not be possible for general processes. For 
example, if individual asset prices follow some local volatility models, see Dupire (1994), then 
there are in general no closed form solutions and the PEA method cannot be applied. 

To price basket options for general asset price processes one may study directly the basket 
value and its associated stochastic processes which may contain stochastic volatilities and/or 
stochastic jump intensities and sizes. Dupire (1994) and Derman and Kani (1998) show that 
any diffusion model with stochastic volatility can be replaced by a local volatility model with- 
out changing the European option price and the marginal distribution of the underlying asset 
price thanks to the uniqueness of the solution to the corresponding pricing equation, a parabolic 
PDE. In fact, Gyongy (1986) discovers the equivalence of a non-Markovian model with a Marko- 
vian model and proves that one-dimensional margins of any ltd process can be matched by a 
one-dimensional Markovian local volatility process, that is, the value of the square of the local 
volatility is equal to the expectation of the square of the stochastic volatility conditional on 
the final stock price being equal to the strike price. Without the uniqueness one cannot claim 
the equivalence of the stochastic and local volatility models but can still derive the same pric- 
ing equation. Bentata and Cont (2009) derive a discontinuous analogue of Gyongy's result for 
semimartingale asset price processes. 

The pricing equation for general asset price processes may contain coefficients expressed in 
terms of some conditional expectations. It is in general a challenging task to compute these 
conditional expectations as there is no closed form solution to the related SDE. One then tries 
to find some good approximations. Antonov et al. (2009) use the Markovian projection onto a 
displaced diffusion and Avellaneda et al. (2002) apply the steepest decent search with Varadhan's 
formula, both methods require to solve some minimization problems. Xu and Zheng (2009) derive 
a closed form approximation to the conditional expectation with a weighted sum of the lower 
bound and the conditional second moment adjustments. 

In this paper we discuss the European basket options pricing for a local volatility jump- 
diffusion model. The main idea is to reduce a multi-dimensional local volatility jump-diffusion 
model problem to a one-dimensional stochastic volatility jump-diffusion model, then to derive 
a forward PIDE for the basket options price with an unknown conditional expectation, or local 
volatility function, and finally to apply the asymptotic expansion method to approximate the 
local volatility function. The main contributions of the paper to the existing literature of the 
basket options pricing are the following: we propose a correlated local volatility jump-diffusion 
model for underlying asset price processes and derive a forward PIDE for general asset price 
processes with stochastic volatilities and stochastic jump compensators, which may be used for 
other applications in pricing and calibration, and we find the approximation of the conditional 
expectation with the asymptotic expansion method. Numerical tests show that the method 
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discussed in the paper, the asymptotic expansion method, performs very well for most cases in 
comparison with the Monte Carlo method and the PEA method discussed in Xu and Zheng 
(2009). 

The paper is organized as follows. Section 2 formulates the basket options pricing problem, 
reviews some pricing results on jump-diffusion asset price models and forward PIDEs. Section 
3 discusses the computation of the conditional expectation and applies the asymptotic expan- 
sion method to approximate the local volatility function. Section 4 elaborates the numerical 
implementation and compares the numerical performance of different methods in pricing basket 
options. Section 5 is the conclusion. The appendix contains the outline of the derivation of a 
forward PIDE for a general stochastic process. 



2 Forward PIDE for Basket options pricing 

Assume a portfolio is composed of n assets and the risk-neutral asset prices Si satisfy the following 
SDEs: 

= (r - \m)dt + Gi{t, S i (t-))dW l (t) + (e Y - l)dN(t) (1) 

where Wi(t) are standard Brownian motions with correlation coefficients pij between Wi(t) and 
Wj(t), (e Y — l)dN(t) is a differential form of a compound Poisson process £3(e y *-l) with N 
being a Poisson process with intensity A and {Yfc} iid normal variables with mean r] and variance 
7 2 , e Yk — 1 is the proportional change of the asset prices at the kth jump of the Poisson process, 
m = Efe^ 1 — 1] = e' 7+ 2 7 — 1, <Ji(t, S) are bounded local volatility functions, r is constant risk-free 
interest rate, and N, W and are independent to each other. The basket value at time t is 
given by 

n 

S(t) = J2wiSi(t), 
i=i 

where Wi are positive constant weights. The basket value follows the SDE 

= (r - Xm)dt + "Wife Si(t-))^^dWi(t) + (e y - l)dN(t). (2) 

Define 

1 n 

V ^ 2 ■= s(tj2 E ^M^S l {t))a^S 3 (t))S i (t)S J {t) Pl] (3) 



then ([21) becomes 

= (r - Xm)dt + V(t)dW(t) + (e Y - l)dN(t) (4) 



S(t- 

with the initial price S(0) = 2^=1 Wj5j(0) and a standard Brownian motion W. Note that V(t) 
is a stochastic volatility which depends on individual asset prices, not just the basket price. We 
have a stochastic volatility jump-diffusion asset price model for the basket option problem. We 
next review some related pricing results for jump-diffusion asset price processes. 

Andersen and Andreasen (2000) model the risk-neutral asset price S by a jump-diffusion 
process 

Jpl = (r(t) _ q{t) _ \(t)m(t))dt + a(t, S(t-))dW(t) + (e y W - l)dN(t) (5) 
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where W(t) is a Brownian motion, N(t) a Poisson process with deterministic intensity \(t), Y(t) 
a sequence of independent random variables with time dependent density <^(-; i), a a bounded 
time and state dependent local volatility function, m a deterministic function given by m(t) = 
E[e Y W — 1], r and q deterministic risk-free interest rate and dividend yield, and N, W, Y 
independent to each other. Andersen and Andreasen (2000) derive a forward PIDE for the 
European call option price C(T,K) at time as a function of maturity T > and exercise price 
K > 0: 

C T (T,K) = -q(T)C(T,K) + (q(T)-r(T) + X(T)m(T))KC K (T,K) 
+ ^a(T,K) 2 K 2 C KK (T,K) 

+ A(T) (J™ C(T,Ke-y)ey((y;T)dy-(l + m(T))C(T,K)\ (6) 

with the initial condition C(0,K) = (5(0) — K) + . Andersen and Andreasen (2000) also discuss 
the stochastic volatility jump-diffusion model and point out that the European call option price 
satisfies the same PIDE (JBJ with the local volatility function a replaced by 

cr(T, K) 2 = B[V(T) 2 \S(T) = K] (7) 

where V(t) is a stochastic volatility process. 

Carr et al. (2004) generalize the asset price model of Andersen and Andreasen (2000) to a 
general local volatility and local Levy type model given by 

= ( r (t) - q{t))dt + a(t, S{t-))dW(t) + f°° (e x - l)[n(dx, dt) - v(dx, dt)] (8) 

where fj,(dx, dt) is a random counting measure and u(x, t) a compensator with a form u(dx, dt) = 
a(S(t),t)k(x)dxdt, here a(S,t) is a deterministic local speed function that measures the speed at 
which the Levy process is running at time t and stock price S, and k(x)dx specifies the arrival 
rate of jumps of size x per unit of time, all other terms are the same as in equation ([5]). Carr et 
al. (2004) derive a PIDE for the European call option price C(T, K) at time as 

C T {T, K) = -q(T)C(T, K) + (q(T) - r(T))KC K (T, K) + ±a(T, K) 2 K 2 C KK (T, K) 

+ a(T, z)C zz (T,z)zij e ( In -jdz (9) 

with the initial condition C(0,K) = (5(0) — K) + , where ip e is the double-exponential tail of the 
Levy measure given by 



^e(y) 



f™(e x — e y )k(x)dx for y > 



I-oo( eV - e x )k(x)dx for y < 0. 



Kindermann et al. (2008) show the existence and uniqueness of the solution to the PIDE §§§ under 
some continuity and uniform positive definiteness conditions. Carr and Wu (2009) generalize the 
local volatility asset price process ([8j) further to a stochastic volatility asset price process with a 
stochastic jump compensator u (dx, dt) = a(t)k(x)dxdt and a(t) being the stochastic instantaneous 
variance. Carr and Wu (2009) use the model and the fast Fourier transform to value stock options 
and credit default swaps in a joint framework. We can show that the European call option price 
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C (T, K ) satisfies the PIDE ([9]) for a general stochastic process a with the local volatility function 
a given by (J7|) and the local speed function a given by 

a(T,z) = E[a(T)\S(T) = z] (10) 

where E[a(T)|5(T)] is the conditional expectation of a(T) given S(T). The derivation of the 
PIDE ([9]) with the local volatility function ([7]) and the local speed function (fTUj) is given in the 
appendix. 

3 Approximation of Local Volatility Function 

For the basket value process the corresponding local volatility function is given by 

1 n 

a(T, Kf = ^2 E WPijEfo(T, Si(T))aj(T, ^(T))|5(T) = K] 

and (ii(T, Si(T)) = (Ji(T, Si(T))Si(T). Piterbarg (2007) uses the Taylor formula to approximate 
<Tj(T, Si(T)) to the first order with respect to Si(T) at point Fj = 5 , j(0)e rT to get 

a i (T,S l (T))^p l + q l (S i (T)-F l ) 

where pi = &i(T, Fj) and % = Jj^d^T, Fj). We use the same first order approximation to get 

<7i(T, Si(T))<jj(T, Sj(T)) « + Pj qi( S i( T ) ~ F i) + PiQjiSj(T) - Fj). 

If we define a(T,K) 2 = a(T,K) 2 K 2 , then 

<t(T, K) 2 w E WiWjPijPiPjil + tpi(T, K) + <pj (T, K)) (11) 

where ^(T,^) = gE[5<(T) - F|5(T) = X]. 

Without loss of generality we may assume r = (otherwise we can work on discounted asset 
price processes) so Fj = S{(0). To obtain analytical approximation to E[£j(T) — S , j(0)|5(T) = K], 
we use the asymptotic expansion approach related to small diffusion and small jump intensity and 
size, see Benhamou et al. (2009). The perturbation and its purpose are different in this paper. 
In Benhamou et al. (2009) the authors expand a parameterized process to the second order and 
apply it directly to price European options. In this paper we use a different parameterized process 
and expand it to the first order to get the analytic tractability and use it to approximate the 
conditional expectation of stochastic variance. In other words, we use the asymptotic expansion to 
find the unknown local volatility function and then use it in the forward PIDE, while Benhamou et 
al. (2009) use a different asymptotic expansion to a process with a known local volatility function 
and then find the options value directly. Assume e € [0, 1] and define 

dSf(t) = -Xm e St(t-)dt + e&i(t, St{t))dWi{t) + Sf(i-)(e ey - l)dN(t) 

with the initial condition Sf(0) = Si(0), where m e = E[e eY - 1] = e^\^ 2 - 1. Note that 

Sj(T) = Si(T). If we define S^kit) = d I e = 0> then the first order asymptotic expansion 
around e = for Sf(T) is 

St(T)^S i>0 (T)+S itl (T)e. (12) 



5 



We can find Si t o(T) and Si t i(T) as follows: Si t o satisfies the equation dSi t o(t) = with the initial 
condition Sifi(0) = S'i(O), therefore, Si t o(t) = Si(0) for all t. S^x satisfies the equation 

dS iA (t) = -\ v Si(0)dt + Si(0))dWi(t) + Si(p)YdN(t) (13) 

with the initial condition Si t \(0) = 0. Here we have used the result Si t o(t) = «Sj(0). Therefore, 

S itl (T) = -XnSi(0)T + / &i(t, Si(0))dWi(t) + 5i(0) V Y h (14) 

■ / ° i=i 

The asset value S{(T) at time T may be approximated by 

Si(T) = S}(T) « Si t o(T) + S it i(T) = 5<(0) + 5i,i(T) 
and the basket value by 

5(T) « 5(0) +£a; i S,i(r) := S C (T). (15) 

i=l 

Conditional on N(T) = k, the variable Si t i(T), written as Sn(T,k), is a normal variable with 
mean (—XT+k)r]Si(0) and variance J of (i, Sj(0))df +&7 2 iSj(0) 2 , and the variable S C (T), written 
as S C (T, k), is also a normal variable with mean 

Mc (&) = (1 - ATr? + *!7/)5(0) (16) 

and variance 



°c(k) 2 = V \( [ &i(t, Si^a^t, Sj{(l))dt) Pij + k 1 2 S l (0)S j (0) . (17) 



ij=i L Jo 
Therefore, 

E[5i(r)-5 f (0)|5(T) = if] « B[S hl (T)\S c (T) = K] 

oo 

= ^P(N{T) = k)E[S i , 1 (T,k)\S c (T,k) = K\. 
Since Si t \(T, k) and S C (T, k) are normal variables, we can find E[Si t \(T, k)\S c (T, k) = K] exactly 



as 



B[S hl (T, k)\S c (T, k) = K]= E[5i ; i(T, k)} + ^-^(K - fi c (k)) 



Cjjk) 
<?c{kf 



where Ci{k) is the covariance of Si t i(T,k) and S c (T,k), given by 

d(k) = Vo;Jpy( / T <7 4 (t,5 4 (0))<T i (t,S i (0))^)+fc 7 2 S i (0) < S,(0) 
.7=1 L Jo 



We obtain <fi{T,K) in dTTJ as 



Pi ^ ^(/c) 2 
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and <t(T, K) 2 in QJ as 

a(T, K) 2 = a{T) + b(T)K - c(T)S(Q) 

where 

n 

a { T ) = E w i w jPijPiPj 

ij=ifc=o ac[ > KVl Pj J 

C ( T ) = E E ^^iT^ vWPijPW f-C^) + -CiW) (1 - ATr? + to/). 

The European basket call option price C(T,K) at time satisfies the PIDE ([6]), i.e., 
C T (T,K) = XmKC K (T,K) + ^a(T,K) 2 K 2 C KK (T,K) 

/oo 
C(T, Ke- y )e y <j) vn 2 (y)dy - A(l + m)C(T, K) (18) 
-oo 

with the initial condition C(0, K) = (S(0) — K) + , where cr(t, S) is a local volatility function given 
by 

Va(t) + 6(t)5-c(t)5(0) 
a(t,S) = . (19) 

and ^n ~a is the density function of a normal variable with mean 77 and variance / y 2 . 

4 Numerical Results 

In this section we do some numerical tests for the European basket call options pricing with the 
underlying asset price processes ([T]). We use three different methods to facilitate the comparison: 
the full Monte Carlo (MC), the asymptotic expansion (AE), and the control variate (CV) method. 

The MC method provides the benchmark results. We use the control variate technique to 
reduce the variance. In (fT5|) the basket value S(T) is approximated by the first order asymptotic 
expansion S C (T) which is used here as a control variate in MC simulation. The basket option 
price with the control variate S C (T) is given by 

00 

E[(S C (T) - K) + ] = J2 P ( N ( T ) = k M(Sc(T,k)-K) + }. (20) 

fc=0 

Since S c (T,k) is a normal variable with mean fj, c (k) and variance a c (k) 2 , see (|16|) and (|17p . a 
trivial calculation shows that 

ekw) - m _ (^f ) + w*> - *)* 

where is the density function of a standard normal variable and $ its cumulative distribution 
function. 

The AE method is to solve the PIDE (|18p with the approximate local volatility function (|19p . 
We find the numerical solution with the log transform of variables and the explicit-implicit finite 
difference method of Cont and Voltchkova (2005). 
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The CV method approximates the basket value S(T) with a tractable variable S C (T) and 
finds a closed form pricing formula (|20p and (I2ip . This approach is essentially in the same spirit 
as that of Benhamou et al. (2009) with the difference that we only expand to the first order while 
Benhamou et al. (2009) to the second order. 

The following data are used in all numerical tests: the number of assets in the basket n = 4, 
the portfolio weights of each asset Wi = 0.25 for i = l,...,n, the correlation coefficients of 
Brownian motions pij = 0.3 for i, j = 1, . . . , n, the initial asset prices S'j(O) = 100 for i = 1, . . . , n, 
the risk free interest rate r = 0, the dividend rate q = 0, the exercise price K = 100. 

Table 1 

Table 1 displays the numerical results of the European basket call option prices with the MC, 
AE, CV methods. The first column is the maturity (T = 1,3), the second and the third the 
coefficients of local volatility functions (c(S') = aS^ -1 with a = 01,0.2,0.5 and (5 = 1,0.8,0.5), 
the fourth the MC results with standard deviations in the brackets, the fifth the AE results with 
relative percentage errors in comparison with the MC results, the sixth the CV results with errors, 
columns four to six correspond to the case of jump intensity A = 0.3, and columns seven to nine 
are similarly defined for A = 1. The last row displays the average standard deviations of the MC 
method and the average errors of the AE and CV methods. For normal variable Y ~ N^,^ 2 ) 
we set T) = —0.08 and 7 = 0.35. Whenever there is a jump event the jump size is relatively small 
(about 2% of the value lost). The choice of 77,7 and intensity A = 0.3 follow those of Benhamou 
et al. (2009) where the authors claim that these parameters are not small, especially for the jump 
intensity A and the jump volatility 7. 

It is clear that the overall performance of the AE method is excellent. All relative errors are 
less than 0.5% except for the four cases when the local volatility function is c(S) = 0.5. This is 
the case corresponding to the high volatility in the Black-Scholes setting and is irrelevant to the 
maturity T and jump intensity A. This is the phenomenon also reported by other researchers. 
The CV method is not satisfactory with average relative error about 7%. We use the Matlab 
to do the computations. When T = 1 we run 30,000 simulations for each case and repeat 10 
times to get the average value, which is used as the Monte Carlo result. We choose the time step 
size 1/512 and state step size 1/1024 for the explicit-implicit finite difference method, it takes 40 
seconds for the AE method and more than 30 minutes for the MC method. When T = 3 we run 
100,000 simulations for each case and repeat 10 times to get the average Monte Carlo result and 
choose the same step sizes as those for T = 1, it takes 2 minutes for the AE method and more 
than 6 hours for the MC method. The AE method is much faster than the MC method while 
the accuracy is reasonable for most cases. 

Table 2 

Table 2 is similar to Table 1 with the only difference that the mean of Y is r/ = —0.3. 
Whenever there is a jump event the jump size is relatively large (about 21% of the value lost). 
The performance of the AE method is very similar to that in Table 1 with the average relative 
error 0.5%, but the performance of the CV method becomes much worse with the average relative 
error 18%. Since the CV method is similar in spirit to the method of Benhamou et al. (2009) 
there is a possibility that large errors may appear when 77 = —0.3 in Benhamou et al. (2009). 

Table 3 
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Table 3 displays the results with three different methods: MC, AE, and the partial ex- 
act approximation (PEA) method suggested in Xu and Zheng (2009) when the local volatil- 
ity function is a(t, S) = 0.2 and the variable Y is a constant — rj. The purpose of the test is 
to see and compare the performance of the AE and PEA methods. The basic data are the 
same as those in Table 1 and 2. We perform numerical tests for three constant jump sizes 
m = e- v - 1 with rj = 0.25,0.125,0.0625, which results in m = -0.2212,-0.1175,-0.0606, 
respectively. The first column is the jump intensity (A = 0.3,1), the second the jump size 
( m = -0.2212, -0.1175, -0.0606), the third maturity (T = 1,3), the fourth the MC results, the 
fifth the PEA results with relative errors compared with the MC results, the last the AE results 
with relative errors. It is clear that both the PEA method and the AE method perform well with 
the relative error less than 1% for all cases, and the former is more accurate than the latter (the 
average relative error 0.1% vs 0.4%). 



Table 4 is similar to Table 3 with the difference that the local volatility function is changed 
to o~(t, S) = 0.5. It is clear that the performance of the PEA method is much better than that of 
the AE method: the former has relative errors less than 1 percent for all cases while the latter 
has relative errors about 2 percent when T = 1 and jumps to about 6 percent when T = 3, 
irrespective to the jump intensities and sizes. We can reasonably say that the PEA method 
is a better approximation method for the European basket call options pricing when the local 
volatility functions are of the Black-Scholes type. However, the AE method is much more flexible 
and can handle general local volatility functions (and stochastic volatilities) and general jump 
variables e Y (*) , the two cases cannot be solved with the PEA method for the time being. 

5 Conclusion 

In this paper we have discussed the European basket options pricing for a local volatility jump- 
diffusion model. We have derived a forward PIDE for the European options price with general 
asset price processes. We have applied the asymptotic expansion method to find the approxi- 
mation of the conditional expectation of the correlated stochastic variance. We have shown the 
excellent numerical performance of the AE method in comparison with the Monte Carlo and 
other methods in most cases. The idea and methodology of the paper opens the way for other 
processes and refinements, for example, we may get better approximation if we asymptotically 
expand to the second order or we may introduce individual jump processes or different jump sizes 
for the common jumps. We are currently working on these problems. 



Table 4 



Appendix 

Outline of the proof of the PIDE |P|) with the local functions ([?]) and U0\) . According to Protter 
(2003), Theorem IV.68, 




+ 




[l[S(t-)<K](e x S(t-) - K)+ + l [s{t -)>Kl(K - e x S(t-))+] fi(dx,dt) 
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where L K is a local time at K of process S. Taking the expectation on both sides, using Fubini's 
theorem and the martingale property, we have 

E[(S(T) - K)+] 

= (S(0)-K)+ + J (r(t)- q (t))E[l m>K] S(t)]dt + ^E[L^} 

[l m <K)(e x S(t) - K)+ + l m>Kl (K - e x S{t))+]a(t)k{x)dx 



+ 



T 



E 



dt. (22) 



We have replaced S(t—) by S(t) due to the time integral taken with respect to the Lebesgue 
measure. Differentiating ()22[) with respect to T yields 



dE[(S(T) - K\ 
dT 



(r(T)-q(T))E[l [S{T)>K] S(T)] + 



1&E[L$] 
2 dT 



+ E 



L(T, K, x, S(T))a(T)k(x)dx 



(23) 



where 

L(T,K,x,S(T)) = [l [s{T) < K] (e x S(T) - K)+ + l [s(T)>K] {K - e x S(T)) + ]. 
Since the European call option price at time with maturity T and exercise price K is given by 

C(T,K) = e~ £ r{t)dt E[{S{T) - K)+]. (24) 

we have (Klebaner (2002)) 



E[lr 



where Fgrr) i s the cumulative distribution function of S(T), and 

dF s(T) {K) 



_ dC(T,K) eJ T rm 



dK 



Jo r ^ dt C KK (T,K). 



(25) 



(26) 



Note that the above equation and derivatives are defined in the sense of distribution. If S(T) 
admits a continuous probability density function then C (T, K ) is twice continuously differentiate 
and (1261) holds in the classical sense. Since 



E[(5(T) - K)+] = E[1 [S(T)>K] S(T)] - KE[1 [S{T]>K] ] 
we can combine ((24"]) with ((23]) to yield 

E[1 [S(T)>K] S(T)] = e£*®*C(T,K) - Ke^^^p-. 
We also clearly have 



8E[(S(T) - K)' 



d 



C(T, K)e& r ^ dt + C(T, K)e% r ^ dt r(T) 



dT dT 
Following the same proof as in Klebaner (2002), Theorem 4, we can show that 

= E[V(TfK 2 \S(T) = K}ef° r ® dt C KK (T,K). 



(27) 
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The equation (j27[) and derivatives are denned in the sense of distribution. Klebaner (2002) proves 
(|27p for continuous semimartingale asset price process, it also works for the case here thanks to 
the property of the local time, that it, 

l-T 



/oo rl 
g(K)L^dK = / g(S(t-))d(S c ) t 
-oo JO 



for all positive bounded functions g, where (S c ) v is the quadratic variation of the continuous part 
of the process S. Everything then proceeds exactly the same. We now estimate the last term in 
([23]) . Using Fubini's theorem and the tower property, also noting (fT0j) and (f26|) . yield 

/oo 
[L(T, K, x, S(T))a(T)\ k(x)dx 
-oo 

E [E[L(T, K, x, S(T))a(T)\S(T)]] k(x)dx 
E [L(T, K, x, S(T))a(T, S(T))] k{x)dx 
L(T, K, x, z)a(T, z)dF s ^ {z)k{x)dx 
L(T, K, x, z)a(T, z)e^ r(t)dt C zz {T, z)dzk{x)dx 
lo r ^ dt J a(T,z)C zz (T,z)(j L{T,K,x,z)k(x)dxjdz 



— oo 
oo 



-oo 

oo />oo 



-oo ^0 

OO fOO 



oo J 



e f T r(t)dt J™ o(T) z)Czz (T, z)zip e (in jj dz 



where tp e is the double-exponential tail of the Levy measure k. The last equality follows exactly 
Carr et al. (2004). Substituting everything into (|23|) and simplifying the expression we then get 
the PIDE © with local functions © and (fT0|) . 
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A 


0.3 


1 


T 


a 


P 


MC (stdev) AE (err%) CV (err%) 


MC (stdev) AE (err%) CV (err%) 


1 


0.1 
0.2 
0.5 


1 


5.91 (0.03) 5.91 (0.0) 6.14 (3.9) 
8.14 (0.02) 8.13 (0.1) 8.31 (2.1) 
15.50 (0.04) 15.18 (2.1) 15.52 (0.1) 


11.86 (0.05) 11.83 (0.3) 12.7 (7.1) 
13.25 (0.06) 13.24 (0.1) 13.85 (4.5) 
18.89 (0.05) 18.60 (1.5) 19.25 (1.9) 


0.1 
0.2 
0.5 


0.8 


4.64 (0.02) 4.64 (0.0) 5.03 (8.4) 
5.47 (0.02) 5.47 (0.0) 5.74 (4.9) 
8.11 (0.03) 8.11 (0.0) 8.29 (2.2) 


11.13 (0.03) 11.16 (0.3) 12.43 (11.7) 
11.60 (0.06) 11.58 (0.2) 12.56 (8.3) 
13.25 (0.04) 13.23 (0.2) 13.84 (4.5) 


0.1 

0.2 
0.5 


0.5 


4.06 (0.02) 4.08 (0.5) 4.81 (18.5) 
4.24 (0.01) 4.25 (0.2) 4.83 (13.9) 
4.85 (0.01) 4.85 (0.0) 5.19 (7.0) 


10.96 (0.03) 10.98 (0.2) 12.40 (13.1) 
11.00 (0.02) 11.01 (0.1) 12.41 (12.8) 
11.24 (0.04) 11.26 (0.2) 12.44 (10.7) 


3 


0.1 
0.2 
0.5 


1 


12.18 (0.02) 12.16 (0.2) 12.75 (4.7) 
15.25 (0.03) 15.14 (0.7) 15.65 (2.6) 
27.23 (0.09) 25.64 (5.8) 27.21 (0.1) 


22.94 (0.11) 22.99 (0.2) 24.36 (6.2) 
24.49 (0.13) 24.45 (0.2) 25.78 (5.3) 
33.03 (0.14) 31.55 (4.5) 34.02 (3.0) 


0.1 
0.2 
0.5 


0.8 


10.69 (0.02) 10.68 (0.1) 11.65 (9.0) 
11.64 (0.02) 11.62 (0.2) 12.29 (5.6) 
15.19 (0.03) 15.14 (0.3) 15.62 (2.8) 


22.45 (0.10) 22.51 (0.3) 23.97 (6.8) 
22.81 (0.10) 22.79 (0.1) 24.19 (6.1) 
24.51 (0.07) 24.48 (0.1) 25.76 (5.1) 


0.1 
0.2 
0.5 


0.5 


10.16 (0.02) 10.16 (0.0) 11.53 (13.5) 
10.29 (0.02) 10.29 (0.0) 11.55 (12.2) 
10.92 (0.02) 10.91 (0.1) 11.77 (7.8) 


22.35 (0.04) 22.43 (0.4) 23.90 (6.9) 
22.43 (0.07) 22.44 (0.0) 23.91 (6.6) 
22.51 (0.10) 22.57 (0.3) 24.01 (6.7) 


Average 


(0.03) (0.6) (6.6) 


(0.07) (0.5) (7.1) 



Table 1: The comparison of European basket call option prices with the Monte Carlo (MC), the 
asymptotic expansion (AE), and the control variate (CV) methods. The asset price processes 
are modelled by SDE ([1]). The table displays results with different maturities T, local volatility 
functions <Ji(t,S) = aS^ -1 , and jump intensities A. The numbers inside brackets in the MC 
columns are the standard deviations and those in the AE and CV columns are the relative 
percentage errors in comparison with the MC results. The data used are: number of assets n = 4, 
weights Wi = 0.25, correlation of Brownian motions pij = 0.3, initial asset prices S{(0) = 100, 
interest rate r = 0, exercise price K = 100, normal variable Y ~ N(ri,^ 2 ) with rj = —0.08 and 
7 = 0.35. 
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A 


0.3 


1 


T 


a 


P 


MC (stdev) AE (err%) CV (err%) 


MC (stdev) AE (err%) CV (err%) 


1 


0.1 
0.2 
0.5 


1 


6.99 (0.02) 7.00 (0.1) 8.4 (20.2) 
8.84 (0.01) 8.84 (0.0) 9.91 (12.1) 
15.89 (0.02) 15.62(1.7) 16.51 (3.9) 


15.23 (0.03) 15.28 (0.3) 18.20 (20.0) 
15.76 (0.05) 15.79 (0.2) 18.62 (18.2) 

20.24 (0.04) 20.02 (1.1) 22.34 (10.4) 


0.1 
0.2 
0.5 


0.8 


6.45 (0.01) 6.45 (0.0) 8.13 (26.1) 
6.72 (0.01) 6.73 (0.2) 8.24 (22.6) 
8.83 (0.02) 8.83 (0.0) 9.89 (12.0) 


15.14 (0.03) 15.17 (0.2) 18.09 (19.5) 
15.20 (0.04) 15.23 (0.2) 18.15 (19.4) 
15.75 (0.02) 15.79 (0.3) 18.62 (18.2) 


0.1 

0.2 
0.5 


0.5 


6.43 (0.01) 6.44 (0.2) 8.12 (26.3) 

6.44 (0.01) 6.44 (0.0) 8.12 (26.1) 
6.49 (0.01) 6.49 (0.0) 8.14 (25.4) 


15.11 (0.03) 15.15 (0.3) 18.07 (19.6) 

15.12 (0.03) 15.16 (0.3) 18.07 (19.5) 
15.15 (0.05) 15.19 (0.3) 18.10 (19.5) 


3 


0.1 
0.2 
0.5 


1 


14.70 (0.02) 14.71 (0.1) 17.46 (18.8) 
16.85 (0.03) 16.79 (0.4) 19.11 (13.4) 
27.99 (0.04) 26.51 (5.3) 29.19 (4.3) 


27.00 (0.07) 27.03 (0.1) 31.92 (18.2) 
28.08 (0.06) 28.04 (0.1) 32.93 (17.3) 
35.31 (0.06) 33.92 (3.9) 39.44 (11.7) 


0.1 
0.2 
0.5 


0.8 


14.27 (0.02) 14.29 (0.1) 17.12 (20.0) 
14.48 (0.02) 14.49 (0.1) 17.30 (19.5) 
16.81 (0.02) 16.80 (0.1) 19.09 (13.6) 


26.64 (0.07) 26.74 (0.4) 31.63 (18.7) 
26.82 (0.07) 26.91 (0.3) 31.79 (18.5) 
28.07 (0.06) 28.09 (0.1) 32.92 (17.3) 


0.1 
0.2 
0.5 


0.5 


14.22 (0.01) 14.23 (0.1) 17.07 (20.0) 

14.23 (0.01) 14.25 (0.1) 17.07 (20.0) 
14.31 (0.02) 14.32 (0.1) 17.15 (20.0) 


26.58 (0.04) 26.68 (0.4) 31.58 (18.8) 
26.62 (0.05) 26.69 (0.3) 31.59 (18.7) 
26.71 (0.06) 26.78 (0.3) 31.66 (18.5) 


Average 


(0.02) (0.5) (18.0) 


(0.05) (0.5) (17.9) 



Table 2: The comparison of European basket call option prices with the Monte Carlo (MC), the 
asymptotic expansion (AE), and the control variate (CV) methods. The asset price processes 
are modelled by SDE ([1]). The table displays results with different maturities T, local volatility 
functions <Ji(t,S) = aS^ -1 , and jump intensities A. The numbers inside brackets in the MC 
columns are the standard deviations and those in the AE and CV columns are the relative 
percentage errors in comparison with the MC results. The data used are: number of assets n = 4, 
weights Wi = 0.25, correlation of Brownian motions pij = 0.3, initial asset prices Si(0) = 100, 
interest rate r = 0, exercise price K = 100, normal variable Y ~ N(t],j 2 ) with r] = —0.3 and 
7 = 0.35. 



14 



A 


771 


T 


MC (stdev) PEA (err%) AE (err%) 


0.3 


-0.2212 


1 


7.35 (0.01) 7.35 (0.0) 7.35 (0.0) 


3 


12.93 (0.01) 12.92 (0.1) 12.85 (0.6) 


-0.1175 


1 


6.08 (0.01) 6.08 (0.0) 6.07 (0.2) 


3 


10.57 (0.01) 10.56 (0.1) 10.49 (0.8) 


-0.0606 


1 


5.66 (0.01) 5.66 (0.0) 5.65 (0.2) 


3 


9.83 (0.01) 9.82 (0.1) 9.74 (0.9) 


1 


-0.2212 


1 


10.78 (0.01) 10.77 (0.1) 10.78 (0.0) 


3 


18.64 (0.01) 18.63 (0.1) 18.57 (0.4) 


-0.1175 


1 


7.28 (0.01) 7.28 (0.0) 7.28 (0.0) 


3 


12.65 (0.01) 12.64 (0.1) 12.58 (0.6) 


-0.0606 


1 


6.02 (0.01) 6.02 (0.0) 6.01 (0.2) 


3 


10.45 (0.01) 10.43 (0.2) 10.37 (0.8) 


Average 


(0.01) (0.1) (0.4) 



Table 3: The comparison of the Monte Carlo (MC), the partial exact approximation (PEA), and 
the asymptotic expansion (AE) methods. The asset price processes are modelled by SDE ([I]). 
The table displays the results with different jump intensities A, jump sizes m, and maturities 
T. The numbers inside brackets in the MC columns are the standard deviations and those in 
the PEA and AE columns are the relative percentage errors in comparison with the MC results. 
The data used are: number of assets n = 4, weights Wi = 0.25, correlation of Brownian motions 
Pij = 0.3, initial asset prices Sj(0) = 100, interest rate r = 0, and exercise price K = 100, local 
volatility function <7j(i, S) = 0.2, and jump variable Y a constant. 
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A 


771 


T 


MC (stdev) PEA (err%) AE (err%) 


0.3 


-0.2212 


1 


14.71 (0.01) 14.66 (0.3) 14.42 (2.0) 


3 


25.69 (0.04) 25.44 (1.0) 24.14 (6.0) 


-0.1175 


1 


14.08 (0.01) 14.03 (0.4) 13.79 (2.1) 


3 


24.61 (0.03) 24.39 (0.9) 23.07 (6.3) 


-0.0606 


1 


13.90 (0.01) 13.85 (0.4) 13.61 (2.1) 


3 


24.32 (0.04) 24.11 (0.9) 22.77 (6.4) 


1 


-0.2212 


1 


16.60 (0.01) 16.55 (0.3) 16.32 (1.7) 


3 


28.80 (0.04) 28.55 (0.9) 27.28 (5.3) 


-0.1175 


1 


14.64 (0.01) 14.59 (0.3) 14.35 (2.0) 


3 


25.52 (0.05) 25.28 (0.9) 24.00 (6.0) 


-0.0606 


1 


14.05 (0.01) 14.00 (0.4) 13.76 (2.1) 


3 


24.55 (0.04) 24.33 (0.9) 23.02 (6.2) 


Average 


(0.03) (0.6) (4.0) 



Table 4: The comparison of the Monte Carlo (MC), the partial exact approximation (PEA), and 
the asymptotic expansion (AE) methods. The asset price processes are modelled by SDE ([I]). 
The table displays the results with different jump intensities A, jump sizes m, and maturities 
T. The numbers inside brackets in the MC columns are the standard deviations and those in 
the PEA and AE columns are the relative percentage errors in comparison with the MC results. 
The data used are: number of assets n = 4, weights Wi = 0.25, correlation of Brownian motions 
Pij = 0.3, initial asset prices Sj(0) = 100, interest rate r = 0, and exercise price K = 100, local 
volatility function <7j(i, S) = 0.5, and jump variable Y a constant. 
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